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SUMMARY 


Load separation is the representation of the load in the test records of geometries 
containing cracks as a multiplication of two separate functions; a crack geometry function and a 
material deformation function. In this paper, load separation is 
of several two-dimensional geometries such as: compact tension geometiy, 
bend geometiy, and center cracked tension geometry and three-dimensional geometries such as 
semi-elliptical surface crack. The role of load separation in the evaluation of the TaZT 
parameter /-Integral and the associated factor ri for two-dimensional geometries is discussed 

sin? rr 6 h alS ° dlSCUSSes the theoretical basis and the procedure for using load separation as a 
simphfied, yet accurate approach for plastic / evaluation in semi-elliptical surface crack which is 
a three-dimensional geometry. The experimental evaluation of J, and partieulaX y for ^e e 
dimensional geometries is veiy challenging. A few approaches have been developed in thiT 
regard and they are either complex or very approximate. The paper also presents the load 
separation as a mean to identify the blunting and crack growth regions in the experimental test 

FractureVechltsts F, ” ally ' SePara “° n “ “ Oology in Elastic-Plastic 


INTRODUCTION 


fnr tv T , he Path independent ^integral was first introduced by Rice (1) as a crack tip parameter 
form can “enaf'*'” 8 ° f ° r "° n - linear elastic materials - ™s line integral 



[ 1 ] 


where F is a contour around the crack as shown in figure 1, W is the strain energy per unit 
volume, T is the tension vector perpendicular to T, « is the displacement in the x-direction and ds 
is an element of the contour. The experimental evaluation of / using the line integral can be 
done by placing strain gauges on the outside surface of the specimen along a contour in a plane 
perpendicular to the crack front. By defining the strain, external traction, and displacement in the 
^-direction along this contour, / can be evaluated. However, this technique is usually inaccurate 


pmmmm 


not FfLIMEP 


703 


\j\ 



and impractical. An equivalent, yet more practical / expression was introduced by Rice (1) as. 


J = 

B 


du 


da 


[ 2 ] 



Figure 1 The /-integral contour. 


where U is the potential energy and can be 
measured as the area under the load-line 
displacement record, a is the crack length, B is 
the specimen thickness and v is the load-line 
displacement. This / expression allowed Landes 

and Begley (2) to develop their multispecimen technique to evaluate /. Their approach requires 
the testing of many identical blunt notched specimens with different crack lengths in order to 
develop the U versus a relationship, from which / can be evaluated. Despite the accuracy and 
reliability of this technique, the high cost of specimens preparation and testing made it 
impractical. A single specimen technique was later developed for special geometries and crack 
ranges where the load can be represented in a particular form as will be discussed next. 

Rice reduced the energy rate interpretation / form given in eq [2], for deeply cracked 
bend specimens (3), into: 


B b 


[3] 


where b is the uncracked ligament. This / expression was based on representing the applied 
bending moment M versus the angle of rotation due to the crack presence 0 as. 

0 = / (M/i 2 ) ^ 

The / expression in eq [3] can be obtained from the test record of a single bend specimen. 
Merkle and Corten (4) also developed a single specimen / expression for compact specimen. 
They represented the load using limit load analysis as: 

P = [o y B(b/2)(2eL)lg(Q p ) [5] 

where o v is the yield strength, 6 p is the plastic rotation due to the crack and a is a function of a 
and b. Using this form, they developed / as: 

j = 2 ( 1 +«) U_ [6] 

1+a 2 Bb 

This J expression can be also evaluated by testing a single compact tension specimen. 

The load expressions in eqs [4] and [5] can be described as separable forms. This means 
that the load is represented as a multiplication of two separate functions: a crack geometry 
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function and a material deformation function and can be written as: 
P or M = G(b or a) . H(Q or v) 

Also, the single specimen J form can be generally written as: 



[ 7 ] 


[ 8 ] 


where the modification factor ri is a function of the crack length and the geometry. Ernst and 
Paris (5) proved that the single specimen J form and ri -factor only exist if the load can be 
represented by a separable form. 


Further studies suggested representing J as the sum of two parts; J el and J pl . The elastic 
part, can be evaluated using Linear Elastic Fracture Mechanics as [K 2 /E'] where K is the stress 
intensity factor and E' is the effective Young's Modulus. The plastic part, J Dh can be written, 
according to the ASTM standard test method (6), as: 


J P i = 


1 pi 


Bb 


[9] 


where A pl is the area under the load versus plastic displacement test record and the ri ,-factor is 
considered as: p 


hp/ ~ 2 for bend specimen and 

h 

h p i ~ 2 +0.522 — for compact specimen [19] 


The latter was the linear approximation of Merkle-Corten expression as shown in eq [6], The 
single specimen J pl form in eq [9] is based on the load separable form given as: 

P = G(b/W)H(v p ,/W ) 


where Wis the specimen width and G (b/W) and H(v p /W) are the geometry and deformation 
functions. 


/-integral was first considered as a stationary crack parameter because it is based on the 
deformation theory of plasticity. However, it was later applied to extended cracks (7) if the crack 
growth is of the order of the non-proportional plastic region at the crack tip. This allowed Ernst 
et al (8) to develop a technique to evaluate a J-R curve for precracked specimen with J updated 
for the crack growth. This technique is also based on the single specimen J form. Thus, load 
separation must also exist during crack growth. 


The available load separation expressions which are based on limit load analysis have 
been mainly developed for two-dimensional bending geometries such as bend and compact 
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tension specimens. This limited the single specimen J form to this geometry only. In this paper, 
load separation will be demonstrated in bending as well as tension geometries, in stationary as 
well as growing crack test records, and in two- as well as three-dimensional geometries. It will 
be also presented as a mean to define crack blunting and growth regions. It will be also used in a 
new r\ pl method and to develop a key curve that can provide predicted load-displacement records. 
By extending load separation to three-dimensional geometries, an equivalent single specimen J 
form can be developed. Very few techniques are currently available for the experimental 
evaluation of J for these complex geometries and they are either impractical or inaccurate. 
Finally, the load separation as a methodology in Elastic-Plastic Fracture Mechanics will be 
presented. 


LOAD SEPARATION IN 
TWO-DIMENSIONAL GEOMETRIES 


Stationary Cracks 


In order to study load separation in stationary crack test records, a separation parameter Sy 
is introduced. This separation parameter represents the ratio of the loads P(aJ and P(cij) in the 
test records of two identical blunt notched specimens with crack lengths a, and cij at the same 
plastic displacement v pl . Thus Sy can be written as: 


S U - 


P(a t ) 


N1 


pi 


[ 12 ] 


If the load is separable, eq [12] can be rewritten as: 


G(bJW) H(y pl /W) G(b./W) 

G{b.lW)H(y pl IW)\ Vpi G(b.!W) 


[13] 


For stationary cracks, Sy should maintain a constant value over all of the plastic region. Figure 2 
shows the test records of ten identical bend specimen made of HY 1 30 steel with different crack 
lengths. The specimens are 0.9 inch thick and 2.0 inch wide and have span to width ratio of 4.0. 
They were originally reported in Ref. (9). When the loads in the different test records were 
divided by that of the specimen with a/W=0.75 at the same plastic displacement, the load ratios 
maintained constant values over all of the plastic region except for a small part at the early 
plastic behavior, see figure 3. Figure 4 shows the test records of ten HY130 center cracked 
tension specimens with different crack lengths and figure 5 shows the separation parameter for 
this set of test records. These specimens (9) were also 0.9 inch thick and 2.0 inch wide. It is 
clear that load separation exists in tension geometry such as the center cracked tension as well as 
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Separation parameter 



Figure 2. The test records of ten single edge notched bend specimens. 


Plastic displacement v p i (in.) 



Figure 3. Load separation in the experimental test records of bend specimens. 


707 




Displacement v (in.) 



Figure 4. The test records of ten center cracked tension specimens. 
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Figure 5. Load separation in the experimental test records of center cracked tension specimens. 
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bending geometry such as the bend specimen. Sharobeam and Landes (10,1 1) studied load 
separation in 12 sets of test records and found that load separation is dominant in all of the 
studied cases. Their study included different work hardening materials, different bending and 
tension geometries and wide range of crack length and geometrical proportions. 


Precracked Specimens 


• u j i- nitially ’ We assume ^at load separation also exists in the precracked test records but 
with different geometry and deformation functions than those for stationary cracks test records. 
Thus, the load m precracked specimen and blunt notched specimen test records can be written as: 

P p = G p (b p /W) H p (y pl /W) for precrackedspecimen and 1 - 14-1 


P b G b (b b /W) H b (v p[ /W) for blunt notched specimen 


where the subscript p denotes precracked specimen while b denotes blunt notched specimen The 
separation parameter can be then written as: 


p 

p 


G p (b p /W)H p (y pl /W) 

p b _ 

V 

G b (b b /W)H b (v pl /W) 


[ 16 ] 


s P b represents the ratio of the load in precracked specimen test record to the load in the test record 
of an identical specimen but with a stationary crack. Since G^bfW) is constant because the 
crack is stationary, S pb can be rewritten as: 

8 * = AG P % lw ^ pb (y pl iw) [17] 

where A is a constant and h pb (v p /W) is the ratio between the two deformation functions at the 
same plastic displacement. If several precracked test records are used with the same blunt 
notched record, one gets S pb for each precracked specimen test record as: 

Sj b = A G p (b p /W) hJ b (v pl /W) [18] 

where i denotes the different precracked specimen test records. If the different S pb 's were plotted 
together versus b/W and collapsed into one record, then S pb would be independent of v t and all 
the records would have the same geometry function G p (b/W). Thus, S pb can be written as: 

s Pb =AG P a> P w) rim 


This would also indicate that the load is separable for precracked test records and can be 
represented as a multiplication of geometry and deformation functions. Figure 6 shows the test 
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records of three precracked specimens (E71, E75, E76) and a blunt notched specimen (E73). The 
specimens are compact tension and made of A533B steel with 1.0 inch thickness, 2.0 inch width 
and 20% side grooving. They were originally used in Ref. (12). The initial crack to width ratios 
for the precracked specimens are 0.608, 0.607 and 0.607 and the final crack to width ratios are 
0.656, 0.71 15, and 0.794 respectively. Figure 7 shows S ph versus b/W for the different 
precracked specimen test records on a log-log graph. All of the records collapsed together into 
one which proves the assumption of load separation. Sharobeam and Landes (13) studied load 
separation in four sets of precracked specimen test records of different materials and with large 
crack growth. They indicated that load separation existed in all of the studied cases. 

In precracked test records, load separation can be also used to define both the blunting 
and growing crack regions. Figure 8 shows both b/W and S ph versus v pl for one of the A533B 
precracked specimens. There are three regions for S pb versus v pl . The first region is very small at 
the early plastic behavior and is usually called the non-separable region. Then there is a second 
region where both S ph and b/W almost maintain constant values. This is the blunting region. In 
the third region, both S pb and b/W decrease indicating crack growth. The blunting and crack 
growth regions can be also identified in figure 7. The blunting region looks like a vertical line 
while the crack growth region looks like an inclined line. 


Plastic r| Development 

Based on eqs [2] and [1 1], i) pl can be written as: 

= b/w §G(bW) 

11 pl Gib/W ) d{b!W) 


[ 20 ] 


The geometry function G(b/W) can be obtained using the separation parameter S tj or S pb . The 
stationary crack separation parameter Sy maintains constant value over all of the plastic region 
for specific a, and values. If Sy is developed for different a, values with respect to a single a, 
value, one can construct the relationship Sy versus a/W (or b/W) as: 


Gib t /W) 

G(bj/W) 


\vary 


b i 


A GibJW) 


[ 21 ] 


where A is a constant equal to the inverse of G(b/W). The Sy expression in eq [21] can be then 
used directly in eq [20] instead of G(bAV) to evaluate T] p/ . Figures 9 and 10 show the variations 
of S, versus b/W on log-log graphs for both the bend and center cracked tension sets of test 
records shown before in figures 2 and 4. The power law function fits well the S if b/W 
relationship in both cases. Thus, Sy can be written as: 

S v = A l ( b/W) m t 22 
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Displacement v (in.) 



Figure 6. The test records of three precracked and one blunt notched compact tension specimens. 



Figure 7. Load separation in the precracked compact tension specimens test records. 
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where A, is a constant and m is the power law exponent. From eqs [20], [21] and [22] n is 
equal to the power law exponent m. For the given set of bend specimens, r\ pl has an average 
value of 1 96 and for the center cracked specimens, it is equal to 0.96. The evaluation of n , 
using the load separation technique is discussed in detail in Ref. (13). P 


For the precracked specimens, the separation parameter S ph is equal to the geometry 
function times a constant as shown in eq [19], Thus by constructing the relationship S Db versus 

. •2" geo i met i ry functlon can be developed. Figure 7 shows the variation of S pb with respect 

to b/ W on a log-log graph for the A533B compact tension set. The relationship can be well 
fitted by a power law function with a power law exponent of 2. 12 which is also the value of n , 
for this set as discussed early. When identical A533B compact tension specimens but with 
stationary cracks used to evaluate r\ pt using the stationary crack separation parameter S u , a value 
ot 2.17 was obtained which is very consistent with the precracked specimens results. 


Thus, r\ pl can be obtained using either test records of a set of identical blunt notched 
specimens with different crack sizes or a precracked specimen test record together with a test 
record of an identical blunt notched specimen. 


The Key Curve 


As the geometry function is developed, 
written as: 


an expression for the normalized load, P N , can be 


P n = 


G(b/W) 


= H(y pl /W) 


[23] 


The normalized load is a function of the plastic displacement only and is independent of the 
crack length. The P N versus v pl relationship is commonly called the key curve. Figures 1 1 and 
12 show the key curves for both the bend and center cracked tension sets. Figure 13 shows the 
key curve for the A533B precracked specimen test records together with three test records of 
identical blunt notched specimens. Our studies showed that the geometry and deformation 
functions are the same for both stationary and growing crack test records of same material and 
geometry. The key curve can be used to construct the load-displacement test record of any 
specimen of same material and geometry with any crack length within the range of the key curve. 
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Pn (MPa) 






v pl /W 

Figure 12. The key curve for the HY130 center cracked tension specimens. 



Vpi/W 

Figure 13. The key curve for the A533 precracked and blunt notched compact tension specimens. 
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LOAD SEPARATION IN 
THREE-DIMENSIONAL GEOMETRIES 


The experimental evaluation of / in three-dimensional geometries is usually very 
challenging. The single specimen / form given in eqs [8] and [9] is not directly applicable to 
three-dimensional geometries because of different reasons. The single specimen / form can be 
used only for single valued / as in two-dimensional geometries. / in three-dimensional 
geometries is not single valued but it varies along the crack front. The single specimen / form 
requires load separation which is not fully studied yet in three-dimensional geometries. The area 
under the test record in the single specimen / form is based on the Load-Line Displacement 
(LLD). In some three-dimensional geometries such as surface cracks, the LLD is insensitive to 
the crack size and the Crack Mouth Opening Displacement (CMOD) is the displacement that is 
commonly measured. The crack in the single specimen / form is defined by one parameter, its 
length, while in three-dimensional geometries, it is defined by two parameters, its depth and 
width. Thus, in order to simplify / evaluation in three-dimensional geometries and develop a 
form equivalent to the single specimen / form in two-dimensional geometries, these issues need 
to be addressed. Load separation in three-dimensional geometries need to be studied. The 
variable / needs to be replaced by an equivalent single-valued /. The two-parameter crack may 
need to be represented by an equivalent single parameter crack. A factor equivalent to r\ pl in the 
two-dimensional analysis needs to be developed. Finally, the relationship between the LLD and 
CMOD needs to be studied. 

We will focus our discussion here on semi-elliptical surface crack as an example of 
three-dimensional geometries, see figure 14a. Both experimental and numerical data will be 
analyzed. The experimental data are for panels made of 2219-T87 aluminum as a base metal 
with tungsten arc- weld seam that is blunt notched with semi-elliptical surface crack. The panels 
are 1.0 inch thick and 4.0 inch wide. These data were originally used by McCabe et al (15) in a 
study on developing /-integral for surface cracks using the equivalent energy approach. The 
numerical data were developed using a finite element model with 2927 nodes and 387 20-nodes 
hybrid brick elements. The mesh in the crack vicinity includes several rings of focused elements 
as shown in figure 14b. The crack tip sides of the first ring elements were collapsed to capture 
the 1/r singularity and the mid-side nodes in the first three rings elements were moved to the 
quarter point to capture the l/r lf2 singularity. The adequacy of the mesh and accuracy of the 
model results were verified by comparing the results with Newman-Raju elastic solutions and 
also some experimental test records, see Ref. (16). The model is used to provide load-LLD 
records, load-CMOD records, and /-integral at different locations on the crack front for different 
combinations of material, crack depth, crack width, and panel size. These numerical data 
together with the experimental data will be used here to address the conditions listed above to 
develop a single specimen / form for surface cracks equivalent to the single specimen / form in 
two-dimensional geometries. 


716 



Load Separation 


Figure 15 shows three experimental test records for different crack sizes and figure 16 
shows the separation results for this set of records. Also, figure 17 shows several numerical 
records for several combinations of crack depth and width and figure 18 shows the separation 
results for this set. a in the numerical test records represents the remote tensile stress and is 
equal to the load per unit area. ft. is clear that load separation also exists in the test records of 
semi elliptical surface cracks and the load can be then written as: 


o = G(a,c)H(v pl ) 


[24] 


where a is the crack depth, c is the crack width and v pl here is the plastic CMOD. 


Effective Crack Length 


311(1 Landes (16,17) studied several expressions for effective crack length and 
concluded that the expression suggested by the R6 method works well for short surface cracks 
1 his expression can be written as: 


7t a 


2 + 


ale 

alt 


[25] 


where a e is the effective crack length and t is the specimen thickness. For deeper cracks this 
expression may need to be modified. 


The Single Specimen Form 


They also developed an energy rate interpretation form for surface crack 


as: 


Jpljav 


dU 


pi 


da r 


pi 


[26] 


where S is the crack front length and J plav is the average J , over the crack front and can be 
written as: 


where 5 is any point along the crack front. From load separation, eq [24], and the energy rate 
interpretation form, eq [26], a single specimen J form can be written as: 


717 








1.1 


01 1.0 H 

u- 
<D 

-O-' 

<u 

E o*9 

o 


o 

Q. 

c 

o 


. 0.8 


D 

CL 

Q) 

CO 


0.7 


0.6 


Semi -elliptical Surface Crack 

B-Q-D-O-O-O-CI-0-E) -B B S B B B B-B-B-D-O-D-B 


, -A A A A A- A*A"^A-ArA-A 


* •►*-*-*- »— i|c 


QPOQD a/t=0.45, a/c=2.0 
4AA4A a/t=0.45, a/c=1.0 
*-*■**-* a/t=0.45, a/c=0.6 

(The dashed lines represent the mean values) 


0.000 


1 — 

0.004 


j — 

0.008 


Vpl/t 


0.012 


0.016 


Figure 16. Load separation in the experimental test records of surface cracks. 
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[28] 


_L_ 3G f 
pl ’ av SG da 


dv 


pi 


where G is the geometry function expressed as a function in a e . The plastic displacement, v pl , in 
eq [28] could be either plastic LLD or plastic CMOD. It can be shown that the ratio between 
both plastic displacements is independent of the amount of plasticity for the same crack size and 
material. Figure 19 shows the ratio between plastic CMOD and plastic LLD at different stages 
of loading represented by amount of plastic displacement, for different crack sizes. The ratios 
maintained constant values over most of the plastic region. Thus, each type of plastic 
displacement can be considered as a constant times the other. Equation [28] can be written in a 
form similar to the single specimen J form in two-dimensional geometries as: 



[29] 


where: 



, geometry, material 


[30] 


The factor C pl is equivalent to the t] pl factor in two-dimensional geometries. Full description of 
C p/ is given in ref. (16). J plav in eq [29] can be obtained from the test record of a single specimen 
In some cases, the maximum value of J on the crack front and not the average value is required. 
Sharobeam and Landes (16) developed a relationship between the maximum J and J plav . They 
found that this relation is almost independent of the material and specimen size and it can be 
represented as a function of a/t and a/c only. 


THE LOAD SEPARATION METHODOLOGY 
IN ELASTIC-PLASTIC FRACTURE MECHANICS 


Figure 20 shows a schematic flow chart of the load separation methodology in Elastic- 
Plastic Fracture Mechanics. By testing a set of blunt notched specimens or at least a single blunt 
notched specimen together with one or more precracked specimens, rip, and the key curve can be 
developed using load separation. Thus the calibration curves J(a,P) can be obtained. The J-R 
curve can be developed using a precracked specimen test record and the value of t|p, obtained 
using load separation. The calibration curves J(a,P) together with the J-R curve provide the full 
elastic-plastic behavior of any structure of same geometry and material as the tested specimen 
even if it has a different crack size. 
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Figure 18. Load separation in the numerical test records of surface cracks. 
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Figure 19. The ratio [plastic CMOD/plastic LLD] versus [plastic CMOD/thickness]. 
















CONCLUSIONS 


Load separation is dominant in both tension and bending two-dimensional geometries It 
can be also extended to three-dimensional geometries such as surface cracks. This Allows the 

exTst Tn 71° r eq “ Single ~^-P-sion for surface cracks. Load section 
exists m both stationary and growing crack test records. It can be used to identify both the 

blunting and growing crack regions. The q p/ factor can be developed using the separation 
p ameters obtained from load separation. Generally, load separation yielded a new simplified 
approach in Elastic-Plastic Fracture. This approach can be applied to any geometry as long as its 
test records are showing load separation. ® 
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